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We discuss how to reveal the massive Goldstone mode, often referred to as the Higgs amplitude 
mode, near the Superfluid-to-Insulator quantum critical point (QCP) in a system of two-dimensional 
ultracold bosonic atoms in optical lattices. The spectral function of the amplitude response is 
obtained by analytic continuation of the kinetic energy correlation function calculated by Monte 
Carlo methods. Our results enable a direct comparison with the recent experiment [M. Endres, 

T. Fukuhara, D. Pekker, M. Cheneau, P. Schaufi, C. Gross, E. Dernier, S. Ruhr, and I. Bloch, 

Nature 487, 454-458 (2012)], and demonstrate a good agreement for temperature shifts induced 
by lattice modulation. Based on our numerical analysis, we formulate the necessary conditions 
in terms of homogeneity, detuning from the QCP and temperature in order to reveal the massive 
Goldstone resonance peak in spectral functions experimentally. We also propose to apply a local 
modulation at the trap center to overcome the inhomogeneous broadening caused by the parabolic 
trap confinement. 

PACS numbers: 05.30.Jp,74.20.De,74.25.nd,75.10.-b 


I. INTRODUCTION 

Collective modes are important for understanding dy¬ 
namic susceptibilities, which include experimentally ob¬ 
served spectral functions and transport properties. The 
situation becomes particularly intriguing in strongly cou¬ 
pled systems, where a simple description in terms of 
weakly interacting excitations is unreliable and perturba¬ 
tion theories fail even qualitatively. Under these condi¬ 
tions, the role of underlying collective modes on dynamic 
susceptibilities becomes increasingly more important but 
is hard to calculate. This is exactly what happens in the 
vicinity of the two-dimensional (2D) Superfluid-to-Mott 
insulator quantum critical point (SF-MI QCP)P. Though 
it is considered to be one of the best studied strongly cou¬ 
pled systems, its quantum critical dynamics is still poorly 
understood, both theoretically and experimentally. 

In superfluids near SF-MI quantum criticality, the ef¬ 


fective field theory in terms of a complex scalar order 
parameter fk features an emergent particle-hole symme¬ 
try and Lorentz invariance, and is expected to have two 
types of collective modes.l^ The first one originates from 
fluctuations of the phase of df and describes a massless 
Bogoliubov-Nambu-Goldstone mode. The second one de¬ 
scribes amplitude fluctuations and is associated with a 
massive Goldstone (MG) mod^, often referred to as a 
Higgs amplitude mode. The fate of the MG-mode in 2D 
is an intriguing and controversial issue because its decay 
into lower-energy gapless modes is found to be strong. 
The mode was argued to become either completely over¬ 
damped (ie., without any resonance type feature in spec¬ 
tral function^ or be detectable as a well-defined reso¬ 
nance peak in certain spectral functions on the super¬ 
fluid side away from (but not on approach to) the critical 

poinPa. 

Recent pr^ress on ultracold atom in optical lattice 
experiment as well as Monte Garlo simulationpfli^. 
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\/N corrections to higher orde i^, a nd non-perturbative 
renormalization group methodJ^^^ have presented solid 
evidence in favor of yet another scenario: A critical reso¬ 
nance in the universal scaling regime. Unlike MG-modes 
in three (and higher) dimensions which become sharper 
when approaching the QCP, the ratio of the width of the 
resonance peak over its mass remains constant in (2-1-1) 
dimensions. 

In order to detect the MG-mode experimentally, the 
scalar spectral function A(a;) (i.e., the correlation func¬ 
tion of I'l'^l) is considered to be the best prob^. On the 
SF side of the transition point, in the scaling limit, A{uj) 
takes the fornP 

Asf(w) (X A^“^/'"<i>sF(^) , (1) 

where A is the Mott insulator gap for the same amount 
of detuning from the QGP, and v = 0.6717 is the cor¬ 
relation length exponent for the U(l) = 0(2) univer¬ 
sality in (2 -I- 1) dimension^^. The universal function 
<i)sF(a:) starts as d)gF(a; —>■ 0) oc and saturates to a 
quasi-plateau with weak w dependence $sF(a^ ^ 1) oc 
2.3-2/i./ ^ ^0.0225 _ intermediate regime between the 
two limits can be constructed numerically, where a well- 
defined resonance peak associated with the critical MG- 
mode is observed at a: = 3. 3(8P. Monte Garlo data 
also suggest that a similar universal resonance, though 
less pronounced, may equally well be seen on the other 
side of the trans ition, a s well as in the quantum criti¬ 
cal normal liquicP^IIlIiSI. These conclusions are yet to 
be confirmed or refuted experimentally. The bottleneck 
of the numerical analysis is analytical continuation of 
data for correlation functions from imaginary to real 
frequencies, A{iujn) —>■ A(w), where iujn = 27rnT with 
n = 0,±1,±2... are Matsubara frequencies. This pro¬ 
cedure is a notorious, ill-posed problem t hat re quires ap¬ 
plication of certain regularization schemej ^^ l ^^ ( and thus 
independent experimental studies are required for final 
understanding. 

The ultracold atom experiment of Ref. [19] aimed at 
detecting the MG-mode in AgF(a;) and confirmed the 
expected softening of the quantum critical spectrum im¬ 
plied by Q but remained inconclusive with regards to 
the existence of a well-defined MG resonance. To obtain 
A{uj), a 2D Bose Hubbard system was gently ’’shaken” by 
modulating the lattice laser intensity (lattice depth) and 
probed by in situ single-site- and single-atom-resolved 
measurements. The observed signal (through temper¬ 
ature increase) featured a broad maximum whose onset 
softened on approach to QGP, in line with the scaling law 
Q. The onset correlates remarkably well with the en¬ 
ergy of the MG-mode, while the ratio of the onset width 
to its frequency was measured to be approximately con¬ 
stant when approaching the critical point. However, a 
resonance-type peak with diminishing width was not de¬ 
tected, which can be interpreted either as evidence for 
the MG-mode being overdamped in the critical state or 
as broadening caused by finite temperature and system 


inhomogeneity (tight confinement) effect^i^. Thus a di¬ 
rect comparison between numerical calculations and ex¬ 
perimental measurements with a common setup is crucial 
to settle the controversy. 

In this paper, we employ an ab initio numerical pro¬ 
cedure based on quantum Monte Carlo simulations and 
numerical analytic continuatiorfi^ to calculate spectral 
functions for ultracold atoms in optical lattices. The fi¬ 
nal result for the temperature increase as a function of 
modulation frequency successfully reproduces the main 
data of Ref. [inifor the experimental setup “as is”; ie, in 
the spirit of the quantum simulation paradigirP^I. The 
consistency between numerical results and experimen¬ 
tal measurements establishes the reliability of both ap¬ 
proaches, and, in particular, validates the analytic con¬ 
tinuation procedure. Moreover, simulations performed 
for various system parameters indicate several improve¬ 
ments/requirements with regards to the experimental 
setup that will help revealing the resonance peak in the 
spectral function. They include (i) the system should be 
effectively homogeneous to avoid inhomogeneous broad¬ 
ening, which can be achieved through confining the lat¬ 
tice depth modulation locally at the parabolic trap cen¬ 
ter; (ii) the detuning from the QCP should be small, 
j/jc < 1.05, where j = J/U is the dimensionless cou¬ 
pling parameter for the Bose Hubbard model introduced 
below [see Eq.([^], and jc is its critical value; and (hi) 
the system’s temperature should be at least as low as the 
Berezinskii-Kosterlitz-Thouless transition point Tc. Our 
results suggest that a direct observation of a well-defined 
resonance peak and understanding the fate of the MG- 
mode experimentally is challenging but not impossible. 

The rest of the paper is organized as follows. We intro¬ 
duce the model in Sec [III and describe the numerical pro¬ 
cedure in Sec. jlllj The comparison between the temper¬ 
ature response from simulations and experimental mea¬ 
surements in a specific setup from Ref. [T9|is presented in 
Sec jlV] We discuss requirements and possible experimen¬ 
tal improvements to reveal the MG resonance in SecjVj 

II. THE MODEL 

Ultracold bosons in optical lattices offer unique possi¬ 
bilities to study the SF-MI quantum phase transition in 
2E| 21 | 22 |. At low enough temperatures the physics of the 
system is restricted to the lowest Bloch band, and can be 
described by the Bose-Hubbard (BH) model 121221^ which 
is parametrized by a hopping amplitude J and an on-site 
interaction energy U, 

Ho = -J '^iblb^+h.c.) + ^'^ni{ni-l)-'^{fj,-Vi)ni, 

(i,i) i i 

( 2 ) 

where (6^) creates (annihilates) a particle on the site 
i, and {i,j) denotes the sum over nearest neighbors on 
the square lattice. In the BH model, the dimensionless 
coupling parameter j = J/U is easily tunable via the 
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lattice depth, and the dimensionality of the system can 
be reduced to 2D by suppressing the hopping in the third 
direction. The total particle number N is controlled by 
the global chemical potential fi. Finally, an ultracold 
atomic gas is trapped by a confining potential Vi, which 
is usually harmonic, Vi = (with m the mass 

of atom, d the lattice spacing, and Ri the distance of site 
i from the trap center measured in units of d). Within 
the local density approximation picture (LDA), ^ — Vi 
plays the role of a local chemical potential. 

Ideally, without the Vi term in Eq. ([^, the system is 
a homogeneous 2D Bose Hubbard model and its phase 
diagram is known with high accuracjl^iHSSI at both zero 
and finite temperature. In the ground state, the system 
undergoes a second order phase transition from the SF 
to the MI when decreasing the ratio j = J/U at fixed 
integer filling factor. At filling factor (n) = 1, the transi¬ 
tion occurs at j~^ = 16.7424(1) and fJ-J/J = 6.21(2) and 
features a QCP with emergent particle-hole symmetry, 
which enlarges the Galilean invariance to Lorentz invari¬ 
ance (the system is actually conformally invariant). The 
SF phase is supposed to have the critical MG-mode ac¬ 
cording to the discussion in the previous section. 

When the Vi term is presented as in current experi¬ 
mental implementationaiSI, due to the inhomogeneous lo¬ 
cal chemical potential, the particle density decreases to 
zero when moving away from the center of the trap. Any 
conclusion regarding the existence of the GM-resonance 
in the homogeneous case cannot be naively applied to the 
realistic experiment, even if the center of the atomic sys¬ 
tem is fine-tuned to be in the vicinity of QGP. A careful 
bottom-up calculation of the scalar spectral function is 
required in order to understand the experimental signal. 


III. SPECTRAL FUNCTION MEASUREMENT: 

THEORY 

In this section, we revisit the generic mathematical 
framework for the measurement of the scalar spectral 
function in ultracold atoms. 

In the BH model, the total kinetic energy K = 
J2{i j) simplest operator with non¬ 

trivial dynamics leading to strong scalar response. Thus, 
we may consider adding an external perturbation term 
Sf{t)K to the Hamiltonian ([^. Within standard lin¬ 
ear response theory, the total kinetic energy response 
is proportional to the external field, and the ratio de¬ 
fines the response function x(w,r) = 6{K{uj))T/Sf{oj) 
where (...)t denotes the thermal average at temperature 
T. The spectral function is defined as the dissipative 
part of the response function, A{uj,T) = 2Imx(w,T), so 
that A{uj,T) is proportional to the energy absorbed by 
the system, which, in turn, determines the temperature 
change of the system. To learn about the spectral func¬ 
tion, one can measure either the total kinetic energy re¬ 
sponse or the temperature change. Though being rather 


indirect, the latter one is the quantity that is measurable 
in the ultracold atom experimeniP^. 

Experimentally, a small uniform modulation 
SVq cos{ujt) of the optical lattice depth Vq is applied 
in t he 2D plane to generate the external perturbation 
ternP^EZl parameter regime where the BH 

model is a valid approximation, the lattice depth in 
units of the recoil energy Er = /2md'^ is much larger 

than unity and controls both parameters J and U: 

J ~ {Vo/ Erf ^ ^ IH 

where the effective dimension D = 2. Substituting Vq 
with Vo + SVo cos{ujt) in J and U, and keeping terms to 
first order in (5Vo; the perturbed BH Hamiltonian reads 

H(t) = Ho + 6g(t)Ho + Sg{t) ^ V^fii + 6f{t)K. ( 3 ) 

i 

Here, the generalized forces 6f{t) = Sfocos{ujt) with 
<5/o = (i - V'^o/Er)^ and Sg{t) = Sgocos{ujt) with 
Sgo = are linear in SVo- Note that the second term 
in Eq. (§ commutes with Ho and yields no contribution 
to the spectral function. Eurthermore, we argue that 
the effect from the confining potential term (third term) 
is also negligible compared to the kinetic energy term 
(fourth term) if one of the following conditions is satis¬ 
fied, i) \5go/dfo\ = 2/\{Ay^Vo/Er-l) \ < 1; ii) the trap is 
large enough and LDA is valid; namely, it is possible to 
decompose the system into independent mesoscopic re¬ 
gions whose sizes are larger than the correlation length 
but are still small enough to be regarded as homoge¬ 
neous systems with the local chemical potential g — Vi. 
This implies, on the one hand, that the total response 
of the system can be approximated by a sum over in¬ 
dependent contributions from mesoscopic regions, while, 
on the other hand, the confining potential term in each 
mesoscopic region is proportional to the particle number 
in this region and thus commutes with the local Ho (i.e., 
it is not dynamic under LDA). The combined effect is 
that the confining potential does not contribute to the 
linear response. The validity of LDA for critical systems 
has been addressed before in Ref. |28| 

For ultracold atoms whose dynamics is dominated by 
the kinetic energy term, the energy dissipation rate is 
proportional to the kinetic energy spectral function A{uj), 

Eico, T) = ^A(a;, T)5fo + P{uj, T). (4) 

Here P{uj,T) is the heating power from other mecha¬ 
nisms (ultracold atoms are always coupled to the photon 
subsystem and are subject to collisions with the back¬ 
ground gas). In the leading approximation, P does not 
depend on the small lattice-depth modulation, and we 
expect P{uj,T) ~ P{T). 

Assuming that the system is quasistatic, i.e., the relax¬ 
ation time is small enough, r ^ 1/w, the thermodynam¬ 
ics can applied at all times and the final temperature 
shift can be deduced from the following self-consistent 
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equation, 

T(w,t)-r(w,o) 


/■* E{u:,T{uj,t')) 
Jo C{T{u;,t')) 


( 5 ) 


where t is the period of the modulation and C(T) is 
the heat capacity. We do not assume here that un¬ 
der the linear response conditions one is allowed to ne¬ 
glect time dependence of temperature on the r.h.s. of 
Eq. (§. This is because over long modulation times 
the temperature change might be substantial. In the ex¬ 
periment^, the initial temperature is chosen to be fre¬ 
quency independent, T{uj,0) = Tjni, and the modulation 
time to be a certain fixed number of modulation cycles 
t = tijjj) = Then the final temperature de¬ 

pendence on frequency, Tgn(a;) = T{u!,t{uj)), is directly 
related to and any sharp resonance structure in 

can be traced back to the spectral function; i.e., 
the temperature response provides a practical probe to 
detect the MG-mode as demonstrated in Ref. m 


IV. MEASURING THE SPECTRAL FUNCTION 
IN SIMULATIONS 

In this section, the experimental setup from Ref. |19| 
is used as a benchmark system for calculation of the 
temperature response from first principles. The pa¬ 
rameters closest to the QCP include the lattice depth 
Vo = lOEr, which gives a dimensionless coupling param¬ 
eter j/jc = 1.2 (or U = 14J), and the particle number 
{N) = 190(36). Combined with the unity filling require¬ 
ment at the trap center, this corresponds to the harmonic 
confinement Vi/J = 0.0915{xi+y/). ^i^The small dimen¬ 
sionless modulation of the lattice depth is SVo/Vq — 0.03, 
which corresponds to the generalized forces \Sfo\ = 0.087 
and |i5(;o| = 0.015. Those parameters define the per¬ 
turbed BH model (|^. 

We argue that the external potential term in Eq. (§ 
is negligible, since both conditions discussed in the pre¬ 
vious section are fulfilled, i) \6go/Sfo\ = 0.17 ^ 1 ; ii) 
the correlation length near the trap center at a typical 
experimental temperature is about one lattice spacin^iS!^ 
so that the LDA also holds in the vicinity of the trap 
center, which dominates the total response. 

In the experiment, the modulation protocol consists 
of two stages. First, ultracold atoms are modulated for 
M = 20 oscillation cycles. Second, the system is held to 
thermalize for some time such that the sum of modula¬ 
tion and hold time is constant at 200 milliseconds for all 
modulation frequencies. During the first stage, both the 
modulation and the heating power P{T) contribute to 
the system’s energy dissipation, while during the second 
stage, only the heating power P{T) contributes to the 
energy increase. The integral over time in Eq. ([^ must 
hence be divided into two segments in order to correctly 
account for both modulation and holding stages. The ad¬ 
vantage of keeping the two-stage time at the same value 
is that the contribution of P{T) is essentially constant 


for all modulation frequencies. Finally, the temperature 
Tfin(w) is determined by slowly ramping the system to 
the atomic limit and measuring the atomic parity in-situ 
with single-site resolution. 

We rely on path-integral quantum Monte Carlo (MC) 
simulations with worm-type updates to calculate the 
scalar spectral functions and the specific heat. Since the 
particle loss in the experiment is negligible, we simulate 
the system in the canonical ensemble. 

In MC simulations, it is straightforward to measure 
the imaginary time correlation function for the kinetic 
energy, x(t) = {K{t)K{0))t — which is related to 

A{uj) via the spectral integral 

p -\-00 7 

Xir) = l ^A^{t,u:-T)A{uj), ( 6 ) 

with the finite-temperature kernel, = 

2(6“'^’’’ -I- We employ the 

same protocol for collecting and analyzing data as in 
Ref. [ini and m More precisely, we collect statistics for 
the correlation function at Matsubara frequencies a;„ 

xiiuJu) = J{\K{iuJn)\'^)T + {K)t, (7) 

and recover x(t) by a Fourier transform. In the path- 
integral representation, x(*w„) has a direct unbiased es¬ 
timator, I where the sum runs over all hop¬ 

ping transitions in a given configuration. Once x('^) is 
obtained from x(*w„) with an accuracy up to 10““^ the 
analytic continuation method described in Ref. is ap¬ 
plied to extract the spectral function A(w). We present 
the analytically continued results at different tempera¬ 
tures in Fig[2 where we see that the curves look quali¬ 
tatively similar at all temperatures in the range between 
0.5 J//cb and 3.33J/fcB; values of A(a;;T) at any temper¬ 
ature in this range can be estimated using linear inter¬ 
polation. All spectral functions vanishing at zero fre¬ 
quency reflects the absence of dissipation in response 
to a static external field. Another way to understand 
this result is through the dissipation-fluctuation theo¬ 
rem, A{uj-,T) = (1 - where S{w,T) is 

the dynamic correlation function of kinetic energy. Zero 
value of A{uj = 0; T) is a natural outcome of a finite 
S'(w = 0;T), see Fig. 6 in Ref. [TUI We also would like 
to point out that the analytical continuation result be¬ 
comes unreliable at very low frequency w <C 1//3 when 
the spectral weight is relatively small. 

The heat capacity for a canonical-ensemble system has 
also been calculated and is shown in Fig. It is seen that 
the heat capacity becomes much smaller in the superfluid 
phase than in the normal phase, which may lead to more 
rapid heating. 

To solve Eq.([^ self-consistently, the initial tempera¬ 
ture Tini and the heating power P(T) are also required. 
However, both quantities were not addressed by the pre¬ 
vious experiment nor is P{T) computable by Monte Carlo 
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FIG. 1. Spectral functions of the kinetic energy per atom at 
different temperatures for the experimental setup of Ref. m 
All curves look qualitatively similar. 


simulations. Thus, we are forced to cousider both quau- 
tities as fitting parameters. In Fig. we show two pos¬ 
sible temperature responses to modulation obtained by 
solving Eq. which both fit the experimental data 

well despite having rather different (but realistic) sets 
of (TinijF). Excellent agreement between the numerical 
and experimental results not only ensures that the an¬ 
alytical continuation procedure (as routi nely app lied on 
the kinetic energy correlation functiorff^^i^E^Ell) is reli¬ 
able, but also validates various assumptions made in the 
first-principle calculation, such as quasi-static thermody¬ 
namics. We also would like to point out that there are 
no fundamental difficulties in experiment to measure Tini 
and P{T), and thus an even more stringent test avoiding 
any fitting can be envisioned in the future. 


C(JI1)IJ peratom 



TU 

FIG. 2. Heat capacity C{T/J) per atom as a func¬ 
tion of temperature T/ J. For the homogeneous system, the 
Berezinskii-Kosterlitz-Thouless transition temperature is at 
fcBlc ~ 1.04J. ^ Notice that below Tc, the heating of the 
system gets boosted due to the smallness of C{T/J) . 


FIG. 3. Temperature response to the lattice-depth mod¬ 
ulation that reveals the spectral function for the ultracold 
atomic system. The vertical axis represents final tempera¬ 
ture. Filled circles connected by a black dashed line are fi¬ 
nal temperatures measured experimentall}^^. The solid lines 
(red and blue) are two predictions based on numerical calcu¬ 
lations, with parameters (feBllni = 0.56J, P = 0.27J/sec) and 
(feBllni = 1.12J, P = 0.45J/sec) respectively. In the calcula¬ 
tions, we assume that the heating power P is independent of 
temperature. 


V. ON THE EXPERIMENTAL OBSERVATION 
OF THE MASSIVE GOLDSTONE PEAK 

Comparing the temperature response in Fig. with 
the spectral function for a homogeneous system with 
j/jc = 1.2 (or C/ = 14 J) (see Fig. 2 of Ref. HO]), we find 
that while the steep onset of the spectral weight corre¬ 
lates remarkably well with the GM-mode energy, the res¬ 
onance structure is lost in the experimental system. As 
mentioned previously, this may occur for two unrelated 
reasons: either because the MG-mode is overdamped, or 
the resonant signal is broadened by finite temperature 
and system inhomogeneity (tight confinement) effectJ^^. 

Our simulations indicate that the second scenario is 
far more likely. Previous work on the homogeneous case 
established that a detuning smaller than j'/jc = 1-05 (or 
U = 16) is required to clearly see the MG-resonance on 
top of the high-frequency continuum. Let us therefore 
take a system with j jjc = 1.05, particle number N = 800 
and unity filling factor at the trap center and perform a 
numerical thought experiment: In order to reduce in¬ 
homogeneity effects, we limit the lattice-depth modula¬ 
tion to a mesoscopic area around the trap center, where 
the confining potential is nearly flat. For simplicity, we 
choose a square area with side length R. In Fig|^ we 
show spectral functions for different values R. Resolving 
the resonance structure hiding in the inhomogeneous sig¬ 
nal is dramatically improved when the modulation area is 
reduced. Though no resonance structure is seen when the 
entire system is modulated {R/d = oo), it emerges when 
the modulation region is reduced to R/d = 16 or R/d = 8 
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at low enough temperature T T,. (where ~ 0.45 J/fee 
is the BKT temperature for a homogeneous model with 
j/jc = 1.05). Converting spectral density to tempera¬ 
ture response does not change this observation qualita¬ 
tively, see Fig. even though the contrast for observing 
the resonance feature is diminishing. This thought ex¬ 
periment demonstrates that by taking care of response 
homogeneity, detuning from the QCP, and temperature, 
the MG-peak can be seen in the kinetic energy spectral 
function using existing technology. 

A{hcolJ)IJ - kJ/J-=3.i3, - k^T/J=2.00 


per atom - kJ/J=l.m - kJ/J=0J0 



FIG. 4. Spectral functions for different modulation areas 
(a square of size R x R whose center coincides with the trap 
center) and temperatures for a system with with N = 800 
atoms and j/jc = 1.05. The MG-resonance emerges at tem¬ 
peratures T ~ J/kB when the modulation is limited to a 
mesoscopic area of linear size R — 16d, where d is the lattice 
spacing. 

Combining the numerical results for the homogeneous 
and inhomogeneous model from Ref. dHI and [T^] and from 
this work, we deduce that the following three conditions 
are to be met in order to reveal the MG-peak: 

First, the modulation area should be restricted to the 
region with nearly constant chemical potential to ensure 
that the temperature response is measured for a homoge¬ 
neous system. To achieve it, a straightforward approach 
would be to replace the harmonic confinement with the 
flat-bottom plus sharp walls potential. This approach, 
however, may lead to problems with controlling system’s 
density and entropy, and , thus, detuning from the QCP. 
An alternative approach is to restrict modulation to a 



FIG. 5. Temperature response for different modulation ar¬ 
eas for the same sets of model parameters as in Fig We 
assume the system’s heating power to be P = 0.2J/sec and 
the initial temperature kBT/J = 0.5. To optimize the con¬ 
trast, the number of modulation cycles is set to be three and 
the sum of modulation and hold time is kept constant at 19 
milliseconds. The best modulation strength are found to be 
SVo/Vo = 0.02, 0.03, 0.06 for R/d = 8,16, oo respectively. 


mesoscopic area around the trap center, as shown is done 
in Fig. One promising experimental implementation 
would be to apply a localized modulation of the scatter¬ 
ing lengtfP^ using a laser beam induced Feshbach reso- 
nance^II. The technique has recent ly be en shown in ul¬ 
tracold atom without optical lattic4221331_ The size of the 
modulation area can be engineered by tuning the size of 
the laser beam (e.g. using a mask). Such a modulation 
would result in a time-dependent on-site interaction U 
in the modulation area. As pointed out in Sec. |III[ by 
subtracting a term proportional to Hq^ the perturbation 
in potential energy can be replaced with the perturba¬ 
tion in kinetic energy, meaning that the MG-mode can 
be studied using the same temperature-response protocol 
as in the current experiment. 

Second, the system has to be close enough to the QGP 
so that a Lorentz invariant action provides an adequate 
description of physics. Our simulations indicate that 
a detuning j'/jc ^ 1-05 is sufficient to reveal the MG- 
resonance, while a smaller detuning j/jc ^ 1.02 is re¬ 
quired to recover the universal spectral function Eq. Q 
including the critical pseudo-plateau at large frequen- 
ciePl. 

Third, the system temperature has to be low enough. 
For a homogeneous system, simulations suggest that the 
resonance peak survives at temperatures as high as the 
BKT transition temperature Tc, but gradually goes away 
at T > 2rE5l. Thus, having initial temperatures below Tc 
is recommended. For example, in the test system with 
j/jc = 1.05, N = 800, and localized modulation size 
R = Sd, which heats from Tc up to 2Tc ( Tc ~ 0.45J/fcB), 
the resonance peak will remain visible in temperature 
response according to Fig. |^C). 
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To conclude, we would like to point out that the quan¬ 
tum critical dynamics in the MI and normal liquid phases 
is also of great interest. Numerical simulations indicate 
the presence of a universal resonance structure in the 
spectral function not only in the SF phase but also in 
phases with un-broken 17(1) symmetry, and at tempera¬ 
tures T 3> 77 (normal quantum critical liquid)?^. The ex¬ 
istence of such universal resonances is unexpected within 
the current weak-coupling theory and their nature re¬ 
quires further study. Verification of this prediction from 
ultracold atom experiments would be crucial to solve this 
puzzle. 
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